Effective separability of typical entangled many-body states 
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We consider two systems of harmonically trapped particles in a typical pure state of the Hilbert 
space defined by given values of the particle numbers and energies of the two gases. Such a state 
is entangled but we show that, for large systems, the resulting correlations between the two gases 
are identical to those of a separable mixture. This result can be generalized to other physical 
systems. We discuss the relation of this effective separability to the well-known existence of quantum 
correlations in any entangled state. We study in detail a small bipartite system and find that its 
correlations are well explained by the large systems results. 
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According to the postulate of equal a priori probabili- 
ties, an isolated quantum system in equilibrium must be 
described by a microcanonical mixed state, i.e., a density 
matrix which commutes with the system Hamiltonian 
and whose only nonzero diagonal elements correspond to 
eigenenergies in a given energy interval. Recently, it has 
been shown that, for many physically relevant properties, 
a pure state description is also possible. An equilibrium 
state is characterized by its energy and possibly other 
thermodynamic variables such as particle number or vol- 
ume. Almost all pure states in the Hilbert space defined 
by these parameters lead to expectation values of inter- 
est identical to those of the corresponding microcanonical 
mixture, provided the considered system is large enough 
[3, 0, HI . Most of the effort has been devoted to deriving 
the canonical ensemble for a system S weakly coupled to 
a large heat bath B from pure states of the composite 
system consisting of S and B @, 0, Q ■ But the pure 
state description of equilibrium is not restricted to the de- 
grees of freedom of a subsystem of a larger system. For 
example, the density profile of a gas in a pure state of 
macroscopically well-defined energy is indistinguishable 
from that of the corresponding microcanonical state Q . 

The effective equivalence of pure states with micro- 
canonical states raises an interesting question in the con- 
text of multipartite systems. Consider two systems A and 
B in a typical state of a Hilbert space characterized by 
some parameters such as the energies of A and B. If the 
above discussed equivalence applies, this entangled state 
cannot be distinguished from the corresponding micro- 
nanical state. But this microcanonical state is separable, 
i.e., is a mixture of product states, showing possibly only 
classical correlations between A and B. On the other 
hand, it has been proved that the entangled character 
of any pure state can be revealed by local measurements 
on systems A and B, irrespective of their nature and 
size @, Q. The o nly pure states which do not violate 
Bell's inequalities [lfj, [llj are product states and hence 
the considered entangled state of A and B must exhibit 
quantum correlations. This seems to be in contradiction 



with a possible equivalence to a separable state. 

In this Letter, we consider two harmonically trapped 
gases of bosons or fermions, in an entangled state charac- 
terized by the energies and particle numbers of the two 
systems. We first show that, for large systems, almost 
all pure states determined by these thermodynamic pa- 
rameters, lead to the same bipartite correlations. These 
correlations are found to be identical to those of a separa- 
ble mixture. For small particle numbers and energies, the 
model we study is simple enough to allow the complete 
determination of the basis spanning the corresponding 
Hilbert space. It is thus possible to evaluate correlations 
between systems A and B for particular entangled states 
drawn from this space. We find that these correlations 
are well described by the expressions obtained for large 
systems, even for as few as ten particles. We discuss in 
detail the relation of our result to the well-known viola- 
tion of Bell's inequalities mentioned above. 

The system A we consider consists of particles confined 
in a harmonic trap and is described by the Hamiltonian 
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where uja is the frequency of the harmonic confining po- 
tential and c Ak creates a particle in the single-particle 
eigenstate tel. Throughout this paper, we use units 
in which H= 1. The cases of bosons and fermions will be 
treated simultaneously in the following. We assume that 
the system A contains a number of particles Na and has 
an energy Ea- The corresponding eigenstates |{nAfe}) of 
Ha satisfy 
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where Ma is an integer, riAk € N for bosons and riAk € 
{0, 1} for fermions. The Hilbert space spanned by these 
states is denoted by Ha and its dimension by Da- The 
system B is described by a Hamiltonian of the form ([I]) 
and is characterized by a particle number Ng and an 
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energy E B = Mb + N B /2. We denote the corresponding 
eigenstates and Hilbert space by |{nsfc}) and Hb- 

For some practical purposes, it is useful to rewrite the 
conditions {2]) as 
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where the positive integers ki obey k%+\ > ki for bosons 
and ki+x > ki for fermions. It is clear from this 
form that Ma can be as small as zero for bosons but 
Ma > Na(Na — l)/2 for fermions. Another interest- 
ing conclusion can be drawn from as follows. A 
fermionic configuration {ki} satisfying with the num- 
bers Na and Ma corresponds to a bosonic configura- 
tion {k'i = ki — i + 1} satisfying Q with the numbers 
N' A = N A and M' A = M A - N A (N A - l)/2. Conse- 
quently, the Hilbert space dimension Da is the same for 
a fermionic system characterized by the numbers Na and 
Ma and for a bosonic system characterized by Na and 
M' A - This dimension increases with Na and Ma- It is 
equal to 1 for Na = 1- For Na = 2, it can be easily 
shown from §5§ that, for bosons, Da — Ma/2 + 1 for 
even M A and D A = M A /2 + 1/2 for odd M A . The di- 
mension in the fermionic case can be obtained using the 
fermion-boson correspondence just discussed. For larger 
Na, Da increases much faster with Ma, see Fig. [2j The 
results shown in this figure are obtained for Na = 5 and 
are well approximated by ln(D^) ~ \J Ma ■ For Ma = 30, 
we find Da = 674. A complete discussion of Da for large 
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Na and Ma can be found in Ref. 

We consider that the bipartite system consisting of A 
and B is in a pure entangled state 

l*> = E ^{n Ak ,n Bk }\{nAk})®\{n B k}) (4) 

{riAk ,n B k} 

which belongs to the product space H = Ha <8> Hb of 
dimension D = DaDb- We are interested in the expec- 
tation values 
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where Ox (X = A or B) is an observable of X with 
eigenvalues between —1 and 1 in a ii^-invariant space 
H' x containing Hx [13]. An example is the number of 
particles of system X between two given positions di- 
vided by the total number Nx 0] • Bell's inequalities are 
written in terms of such observables 

BSElHEJ]- we 

now show that, in the limit of large D, the expectation 
value ([5]) is the same for almost all states |^) EH. To do 
so, we use the normalized uniform measure on the unit 
sphere in H 
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FIG. 1: Correlation {O^O^) as a function of q for bosons, 
N A = N B = 5 and M = 20 (squares), 25 (circles) and 30 
(diamonds). For each M, two typical states Q are shown. 
The lines correspond to the Hilbert space average (0. The 
dashed line is this average for fermions, Na ~ Nb = 5 and 
M = 30. This curve is flat for q > 25 as n k = for k > 25 
in this case. In general, the fermionic curve for a given M is 
close to the bosonic curve for M' = M — 10 for < q < M' , 
and is flat for M - 5 < q < M. 



where {n a } stands for {riAk,nBk}- Using the expression 
ir p R 2p /pl for the volume of a 2p-dimensional sphere of 
radius R, we find the Hilbert space average 



(OaO b ) = / d 2U *Li(0 A B ) = (Oa)e a (Ob) 



B)E b 



where d 2D ^> = Y[i n \ rfRe*{„ a }dIm*{ nci }, 

(Oa)e a = ^ E {{nAk}\0 A \{nAk}) 



(7) 



(8) 



{riAk} 



and (Ob)e b is given by a similar expression. For the 
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Hilbert space variance a 2 = (OaOb) 2 — (OaOb) , we 
obtain, by an analogous calculation, 

1 
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The above sums run only over the configurations satis- 
fying An upperbound to the variance a 2 is thus 
obtained by replacing one of these sums by a sum over 
the set of states |{n a }) spanning H' A ® H' B . Doing so, 
we find a 2 < D~ l {0\) e a {0 2 b ) e b < D^ 1 . In conclusion, 
in the large D limit, (OaOb) is given by the product 
{Oa}e a (Ob)e b for almost all states |^) e H. In other 
words, for given observables Oa and Ob, the correlation 
(OaOb) for a typical entangled state (UJ is identical to 
that for the separable state D^ 1 J2{n a } \{ n a})({ n a}\- 

This result seems to be in contradiction with the fact 
that the entangled nature of the state |^) can be re- 
vealed by performing local measurements on A and B. 
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FIG. 2: Variance a 2 as a function of M for bosons, Na = 
Nb = 5, (Oa, Ob) = (0? i Os ) (solid line and open squares) 
and (O^ ' , ) (dashed line and filled circles). We remark 
that a 2 vanishes for M < 4 in the first case and for M < 
1 in the second case. The lines correspond to the Hilbert 
space variance ([9}. The symbols are obtained by averaging 
((OaOb) - (Oa)e a {Ob)e b ) 2 over 10 3 typical states ©. For 
M > 10, the variance for (Oa,O b ) = (O^ ,0$) where < 
p,q < M lies between the two curves shown in the figure. The 
dotted line is D' 1 = D~ A 2 . 

It has been shown that, for any given entangled state of 
an arbitrary bipartite system (A, B), there exist four ob- 
servables Oa, O a , Ob and 0' B with eigenvalues ±1 such 
that > 2 where = O A (0 B + 0' B ) + 0' A (0 B - 0' B ) 
[H, On the other hand, the above argumentation can 
be applied to F and one finds, in the large D limit, 
(F) = (O a ) Ea (Ob+0' b ) Eb + (0' a ) Ea (Ob-0' b ) Eb which 
is clearly between —2 and 2. In short, for any given ob- 
servables Oa, 0' a , Ob and 0' B , \(F)\ < 2 for almost all 
entangled states ((4]) but, for each of these states, there 
exist purpose-built observables for which this inequality 
is violated. For large many-body systems, these very 
special observables might be difficult to measure as the 
number of measurements that can be made in practice 
is limited. For example, measurements of single-particle 
observables ^2 k k , ^kk , c i xk c xk , such as particle number 
densities, can be achieved, whereas n-particle observables 
with n of the order of Nx, are practically inaccessible. 
Another point is worth mentioning here. In the following, 
we consider an observable F such that (F) — > 2 _ in the 

large D limit. In such a case, although (F) 2 — (F) — > 0, 
the proportion of states (|4|) such that (F) > 2, does not 
necessarily vanish in this limit and the above conclusion 
may be not strictly correct. However, the proportion of 
states |\E') such that (F) > 2 + e where e is a positive 
number, as small as we please, obviously vanishes. 

The effective equivalence of a typical state (j4]) with a 
separable state can be generalized to other systems. First 
of all, for two systems A and B described by the Hamil- 
tonian Jl}, other Hilbert spaces Hx can be considered. 



FIG. 3: Distribution P ((OaOb)) as a function of a = 
{(OaOb) - (Oa)e a (O b )e b )/o for bosons, N A = N B = 5, 
(O a ,O b ) = (C^ 0) ,C4 0) ) (open squares) and (0^ M) ,0^ M) ) 
(filled circles), and M = 10 and 20. The dotted lines cor- 
respond to the normal distribution of variance a 2 . 



The above derivation remains valid, for example, if the 
condition ^ is replaced by Yli=i < Ma or, in other 
words, for states |\&) given by @ with the sum running 
over all the eigenstates ^n^fc}) and \{riBk}} correspond- 
ing to eigenenergies lower than Ea and E E , respectively. 
For other systems, such states show effective separability 
in a proper thermodynamic limit. It can be seen as fol- 
lows. In this limit, the average © and Sa — ks Ih(Da) 
where k E is the Boltzmann constant, are, respectively, 
the usual microcanonical average and entropy [l5| . As 
is well known, the entropy Sa is extensiveand hence the 
variance ([9]) vanishes exponentially with the system size. 

We now study in detail a small bipartite system. We 
consider two systems described by the Hamiltonian ([T]) 
with N A = N B = 5 and M A = M B = M between and 
30 for bosons and between N A (N A - l)/2 = 10 and 40 
for fermions. The corresponding eigenstates are deter- 
mined using As observables Ox, we use the opera- 
tors defined by 0^\{nxk}) = ±\{ n xk}} with the upper 
sign if nx q = where q is a given positive integer, and 
the lower if nx q ^ 0. For these observables, the first 
term of ([9]) simplifies to 1/(D + 1). We remark that 

Ojp = 1 — 2c Xq c Xq for fermions. To draw a normalised 
state Q from the uniform distribution ©, we gener- 
ate D random complex numbers ${ na } with standard 
normal distribution and then compute the components 

= $ {« Q }/E{„ Q } l $ {n Q }| 2 M- For each state ' we 
evaluate correlations ([S]). Results obtained in this way 
are shown in Fig. [T] The agreement with the micro- 
canonical average ([7]) is excellent. As discussed above, 
the dispersion of (OaOb) around this mean value de- 
creases with increasing M, see Fig. [21 Distributions of 
correlations (0^0 { p) are shown in Fig. [3] They are 
constructed by evaluating such expectation values for 10 7 
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evaluated. Then a residue calculation gives 

D e -x(V2+l) 



FIG. 4: Distribution P{(F)) as a function of D((F) - 2) for 
bosons, Na ~ Nb = 5 and M — 2 (squares), 3 (circles) and 
10 (diamonds). The line is the large D expression 



typical states It can be shown that they converge to 
normal distributions as M is increased. 

To discuss Bell's inequalities, let us consider two par- 
ticular states Kn^}) satisfying and define observ- 
ables a z A and a A which act on the states |{TMk}) ^ 
\{n\ k }) like the identity operator, and are represented 
by the Pauli matrices a z and a x , respectively, in the 
basis K^Afc})} °f the corresponding subspace. 

For system B, we define the observables a B and a B 
which coincide with the unit operator except on the sub- 
space spanned by two states K^fc}) m which a' B — 
2 ~ 1/2 (°b - °%) and CT s = 2_1/2 (cr| +a%) where a z B and 
a B are the analogues of a z A and a A , respectively. We are 
interested in the expectation value 



(F) = (a%(a B +a B )+a%(a\ 



,)) (10) 



in a typical state It can be written in terms of the 
four components ^{ na } corresponding to { n Bk})- 
To simplify the following expressions, we denote these 
components by Sfr±±. We first observe that, for states 

|*) such that * ++ = * = and = 

(F) = 2V2r] + 2(1 - rj) where r) = 2|*+_| 2 e [0,1], 
and hence varies between 2 and 2y/2 which is the max- 
imum possible value for such an expectation value [3]. 
So, the CHSH inequality pj, \{F}\ < 2, is violated by 
some states ([JJ. For the Hilbert space average, we find 
(F) =2 — 8/D < 2. It is thus interesting to study the 
distribution P((F)) resulting from the measure ((6|). Fig- 
ure [4] shows such distributions constructed by evaluating 
(F) for 10 s states \^}. We remark that P is independent 
of the choice of the states it depends only on 

the dimension D. A large D expression for P can be de- 
rived as follows. For large D, the distribution of the four 
components ^±± is essentially Gaussian. Consequently, 
using 5(x) = J dk exp(ifcx)/27r, P can be written as the 
Fourier transform of a Gaussian integral which is readily 
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+ 2e x (x-2) ) e(-x) 
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where x = D((F) — 2)/2. This expression agrees very 
well with the results obtained for M as small as 10, see 
FigfJ] From it, we infer that the proportion of states ((3]) 
such that (F) > 2, is (40 + 28V2)" 1 ~ 0,013. 

In summary, we have studied two harmonically 
trapped gases in an entangled pure state characterized by 
the particle numbers and energies of the two systems. For 
almost all such states, the correlations between the two 
systems are identical to those of a separable mixed state. 
We have proved this for large systems by evaluating the 
Hilbert space average and variance of such a correlation, 
and have shown how this proof can be applied to other 
physical systems. We have also studied in detail a small 
bipartite system and found that its correlations are well 
explained by the large systems results. To discuss the 
seeming inconsistency between the effective separability 
found here and the unavoidable violation of Bell's in- 
equalities, we have considered an observable which leads 
to a maximal violation of the CHSH inequality [ll| for 
some of the considered states. The proportion of states 
which do not satisfy this inequality remains finite in the 
limit of large systems but the inequality violation is less 
and less pronounced as the system size is increased. 
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